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We theoretically study the symmetry properties of the single-band Hubbard model with general 
spin-orbit coupling (SOC) on the Kagome lattice. We show that the global U(l) spin-rotational 
symmetry is present in the Hubbard Hamiltonian owing to the inversion symmetry centered at sites. 
The corresponding spin Hamiltonian has, therefore, the SO(2) spin-rotational symmetry, which can 
be captured by including SOC non-perturbatively. The exact classical groundstates, which we obtain 
for arbitrary SOC, are governed by the SU(2) fluxes associated with SOC threading the constituent 
triangles. The groundstates break the SO(2) symmetry, and the associated Berezinsky-Kosterlitz- 
Thouless transition temperature is determined by the SU(2) fluxes through the triangles, which we 
confirm by finite temperature classical Monte Carlo simulation. 

PACS numbers: 71.10.Fd, 71.70.Ej, 74.62.-c 


I. INTRODUCTION 

The Hubbard modeP- has been acknowledged as 
a paradigm of strongly correlated electron systems. 
Despite its simplicity, the Hubbard model can ex¬ 
hibit various phenomena from antiferromagnetisnP 
to metal-insulator transition 3 and high-temperature 
superconductivity 4 . The Hubbard model without spin- 
orbit coupling (SOC) has nontrivial symmetry properties 
besides the apparent global SU(2) spin-rotational sym¬ 
metry, which has guided us to the exploration of rich 
phases and excitation^®. Yang and Zhang^, for ex¬ 
ample, revealed the SU(2) pseudo spin symmetry, which 
includes the U(l) phase symmetry as a subgroup, and 
predicted the massive collective modes in any phase- 
symmetry-breaking superconductivity. 

In spin systems, SOC gives rise to frustration on spin 
interactions and reduces the symmetry in general. 1 -^ The 
spin-dependent hopping, the manifestation of SOC in the 
kinetic terms of the Hubbard model, can be described by 
an SU(2) gauge field. n ®3ln open-ended one-dimensional 
chains, the SU(2) field can be gauged away by a string of 
gauge transformation^ 11 , wherein the global SU(2) spin- 
rotational symmetry is intact. In rings, however, the 
SU(2) field creates a nonvanishing flux in general, which 
makes the system frustrated and reduces the symmetry 
down to U(1)P SU(2) symmetry is rec overed only when 
the enclosed SU(2) flux vanishes ! 15 ! 16 ! Two-dimensional 
lattices are composed of interconnected loops, each of 
which embraces the flux. General SOC breaks the contin¬ 
uous symmetry^ 7 and engenders a long-ran ged magnetic 
order escaping the Mermin-Wagner theorenpD. 

The geometry of the lattice is another source of frus¬ 
tration. The Kagome lattice, a two-dimensional lattice of 
corner-sharing triangles, is a prototypical example that 
brings geometric frustration to antiferromagnetic materi¬ 
als exemplified by herbertsmithite ZnCii 3 (OH) 6 Cl 2 . The 
Kagome lattice Hubbard model without SOC has been 


extensively studied in metal-insulator transitions 4 ®! 
and the van Hove fillingpftiH. The corresponding spin 
Hamiltonian of the Hubbard model in the large-t/ limit 
at half filling has been studied in search of exotic phases 
on the Kagome lattice, such as spin liquidPJ The phys¬ 
ical effects of SOC in the spin Hamiltonian have been 
studied by including its leading order contribution to the 
Hamiltonian known as the Dzyaloshinskii-Moriya (DM) 
interactio n 10 ! 26 !, which has been known to induce a long- 
ranged magnetic orderP 27 ®3 

In this paper, we show that the global U(l) spin- 
rotational symmetry is present in the single-band SOC 
Hubbard Hamiltonian on the Kagome lattice owing to the 
inversion symmetry centered at sites 28 . The correspond¬ 
ing spin Hamiltonian has, therefore, the SO(2) spin- 
rotational symmetry, which can be captured by includ¬ 
ing SOC non-perturbatively. The exact classical ground- 
states, which we obtain for arbitrary SOC, are gov¬ 
erned by the SU(2) fluxes associated with SOC thread¬ 
ing the constituent triangles. The groundstates break 
the continuous symmetry, and the associated Berezinsky- 
Kosterlitz-Thouless (BKT) transition temperature is de¬ 
termined by the SU(2) fluxes through the triangles, which 
we confirm by finite temperature classical Monte Carlo 
simulation.^ 


II. SUMMARY OF MAIN RESULTS 

We study the the single-band Hubbard model to de¬ 
scribe SOC electron systems on the lattice®! 

Hi = ~t c)U jk c k + U^h^hji, ( 1 ) 

(j,k) 1 

where (j, k) represents the nearest neighbors j and k , 
= (c 4 |,c 4 ^) and Cj = (cj-j-, cyj.) T are the electron cre¬ 
ation and annihilation operators, and hj a = cj Q c ?a is the 
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FIG. 1. (Color online) (a) The site labeling scheme and the 
direction of the links, j —» k, of the DM vectors d jk- (b), (c) 
The DM vector d jk of the gauge fields Ujk = exp (—idjk ■ er/2) 
on the link j —> k, directed counterclockwise around triangles 
(b) in the original Hubbard Hamiltonian Hi and (c) in the 
gauge-transformed Hamiltonian Ho, (for n = z), respectively. 


electron number operator with the spin a. Here, t is a 
real hopping magnitude 3 -^ 


U jk = exp(-idj k ■ a/2) (2) 

describes the effect of SOC, which rotates spin of an elec¬ 
tron while hoppinj|- 37 i djfc is the direction of the DM vec¬ 
tor; er is the vector of Pauli matrices; U is the magni¬ 
tude of the on-site Coulomb repulsion. Hermicity of the 
Hamiltonian requires Ujk = U\j , and thus = —djy. 

The DM vectors d jlc are physical, but can be consid¬ 
ered as a particular realization of the SU(2) gauge field 
in the lattice gauge theorj®], which provides a suitable 
language to study the symmetry of the Hubbard Hamil¬ 
tonian. The local SU(2) gauge transformation, £j i-T VjCj 
and Ujk >->■ VjUjkV^, corresponds to the rotation of local 
spin axes. 

The symmetry of the Hubbard Hamiltonian with SOC 
is closely related to the SU(2) flux vector $ enclosed by 
loops on the lattice E2I2U it i s defined by 

exp(-T$ • er/2) = \\_U jk (3) 

j-*k 

for each loop, where j —► k means that sites are traversed 
counterclockwise around the loop as shown in Fig. 0a). 
In the absence of SOC, the fluxes vanish, which results in 
the invariance of the Hamiltonian under the global SU(2) 
spin rotation, Cj i V6j. A finite SOC causes nontrivial 
fluxes through loops, which would reduce the symmetry 
from the continuous SU(2) to the discrete Z 2 , Cj 1 —> — Cj , 
generally. 

The continuous symmetry, however, can persist even 
when SOC is present, if the SU(2) fluxes through the 
loops meet a certain condition. Our main discovery is 
the U(l) spin-rotational symmetry of the Kagorne lattice 
Hubbard model with SOC. In the Kagorne lattice, the 
inversion symmetry centered at sites 28 demands match¬ 
ing of the DM vectors of adjacent triangles. As a re¬ 
sult all the triangles share the same SU(2) flux structure, 
exp(—■ er/2) = U 12 U 23 U 31 , with the site labeling as 
shown in Fig. 0a) and0b), where = <f>n is the SU(2) 


FIG. 2. (Color online) (a) One of the exact groundstates of 
the spin Hamiltonian Hf Q for the angle 0 < (j> < 2-7t/3 
and n = z. Any global spin rotation of the state about the 
axis n yields also a groundstate. (b) BKT transition tem¬ 
perature of _ff| as a function of the angle (j>. The transition 
temperature Tkt is modulated by the flux through triangles, 
<I> = 30. Dots: finite temperature Monte Carlo simulation re¬ 
sults. Line: theoretical result Tkt (161 with the renormalized 
spin length S = 0.69. 


flux vector threading the triangles. For an isolated trian¬ 
gle, the Hubbard Hamiltonian with arbitrary SOC pos¬ 
sesses the global U(l) spin-rotational symmetry, and it 
can be revealed by the local SU(2) gauge transforma¬ 
tions that are determined by the SU(2) flux through the 
trianglef 11 Sharing the same SU(2) flux between the ad¬ 
jacent triangles in the Kagorne lattice extends the U(l) 
symmetry of an isolated triangle to the entire lattice, 
which becomes visible in the gauge-transformed Hamil¬ 
tonian, 

H 2 = -t £ [cje-^^c'k + H.c.] + UJ2 > ( 4 ) 

j->k j 

characterized by the single DM vector 0n [Fig. [ljc)], 
where c' = UjC 3 is the new electron operator and Uj de¬ 
scribes the local gauge transformation that is governed 
by the SU(2) flux. The angle 0 is uniquely determined up 
to 27 t/ 3 by the SU(2) flux through the triangle, 30 = *I> 
mod 27 t. 

The continuous symmetry also manifests itself in the 
corresponding spin model, 

Hi = J S, • R(d lk ) ■ S fe , (5) 

that is derived from the larg e-U limit of the Hubbard 
Hamiltonian Hi (JT) at half filling. Here J = 8 t 2 /U sets 
the energy scale of spin interactions; R{djk) is the SO(3) 
matrix of a rotation about the axis with the angle 
<pjk = |djfc|. The summand can be split into three terms 
with the aid of Rodrigues’ rotation formula: 

TSj • R^djk') * S k —J COS 0j k Sj • S k 

J sin <j)jk djfc * Sj x Sfc 
+ cos0j/c) (Sj ■ dj/ c )(S/ c * djfc). 

The first term is the antiferromagnetic Heisenberg inter¬ 
action; the second term is the DM interaction, which is 
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antisymmetric in exchanging two spins; the third term is 
the symmetric anisotropic interaction which always ac¬ 
companies the DM interaction in insulators! 11 * 12 * 39 ! Ap¬ 
plication to the spin Hamiltonian H[ of the SO(3) equiv¬ 
alent of the SU(2) gauge transformation connecting H\ 
and H 2 yields the new spin Hamiltonian, 

= ,/ S' • R(<j>n) ■ S' fc , (6) 

j-*k 

which has the transparent global SO(2) spin-rotational 
symmetry that is possible to obtain only if we treat the 
SOC non-perturbatively. 

We shall provide the exact classical groundstates of the 
spin Hamiltonians i?| later, from which the ground- 
states of H[ ([5]) can be obtained by reversing the gauge 
transformation. The groundstates of an isolated triangle 
is obtained by the exact diagonalization of the Hamil¬ 
tonian. Inversion symmetry centered at sites makes the 
groundstates of adjacent triangles compatible, which al¬ 
lows us to minimize the spin energy on the entire Kagome 
lattice. These groundstates break the SO(2) symmetry 
of the spin Hamiltonian. We confirm that the associated 
BKT transition occurs at a finite temperature Tkt mod¬ 
ulated by the flux 3 <j> through the triangle with the aid 
of finite temperature Monte Carlo simulation as shown 
in Fig.gb). 

III. SYMMETRY OF THE SOC HUBBARD 
MODEL 

A. Symmetry on the ID lattice 

For an open-ended one-dimensional chain 
{ci,C 2 ,--- , cjv}, one can keep the first electron ci 
unchanged, and perform a string of successive SU(2) 
gauge transformations on the rest electrons by 

c'j = f II I ^ , (7) 

\i<k<j ) 

which transforms the original Hubbard Hamiltonian 
Hi (0 to 

H = + H.c.] + UJ2 n'jtn'jl > ( 8 ) 

3 3 

where the global SU(2) spin-rotational symmetry is 

evident.™# 

For a one-dinrensional ring, the SU(2) symmetry is not 
present generally because of nontrivial SU(2) flux pen¬ 
etrating the ring. The Wilson line for the ring is the 
product of the link gauge fields, which is given by 

e -i *.„/2 = Y[ u jj+1 (9) 

l<j<N+l 


with the periodic boundary condition cy+i = Ci as¬ 
sumed, where $ = d>n is the SU(2) flux vector. The 
flux magnitude $ is uniquely defined up to 2 ti® 1 , and the 
global SU(2) symmetry is restored only when the flux 
vanishes, $ = 0 (mod 2 tt)P^ The flux vector $ can be 
evenly distributed to all links by the local SU(2) gauge 
transformation, which is given by 

c' = S:i-' )^/ 2 N | Y[ u Kk+ i ) dj . ( 10 ) 

\ 1 <k<3 ) 

The resultant Hamiltonian is H 2 0 with <j> = &/N, 
which is invariant under global U(l) spin rotation about 
the axis n, c' 1 —> e _ien cr / 2 c' for any angle 6 W^ 

B. Symmetry on the Kagome lattice 

We show that the Hubbard Hamiltonian on the 
Kagome lattice is invariant under the global U(l) spin 
rotation, which is protected by the inversion symme¬ 
try centered at sites that are respected in many ma¬ 
terials such as herbertsmithite ZnCu 3 ( OH) 6 Cl 2 and Fe 
jarosite compound KFe 3 ( 804 ) 2 (OH) 6 P — *^* The SU(2) 
link gauge field Uj k = exp(— idj k -(r/ 2 ) can be attributed 
to the electro static potential U(r) induced by surround¬ 
ing molecules ! 43 * 44 * The DM vector d j k oc VF x (r^ — r^) 
is invariant under the inversion r 1 —► —r centered at sites 
provided that V (r) is even under the inversion. For ex¬ 
ample, in Fig. [lja) and (b), under the inversion centered 
at the site 1 , the site 2 (operator c 2 ) maps to the site 
2 ' (operator 62 '), which transforms the associated kinetic 
term: 

— tc\ exp(—idi 2 • <x)c 2 -A —tc.\ exp(—zdi 2 • <x)c 2 ' • ( 11 ) 

The invariance of the Hamiltonian under the transfor¬ 
mation requires that the right side of the equation is 
equivalent to —tcj exp(—zd 12 ' • <t)c 2 /, an d, thus, dictates 
di2 = di2'. 

Each site of the Kagome lattice can be labeled by three 
numbers 1,2, or 3, as illustrated in Fig. 0 a). Once 
the link gauge fields f 7 i 2 ,f 7 2 3 , and U 3 i of an arbitrarily- 
chosen triangle (e.g., a shaded one in Fig. [l]) are fixed, 
all the other link gauge fields on the Kagome lattice are 
determined by the inversion symmetry. Since all the tri¬ 
angles have the same SU(2) link gauge fields, the gauge 
transformations ( | 10 [ ) for neighboring triangles are com¬ 
patible. Specifically, the gauge transformation 

{ Cj , if j is labelled by 1 

e i*..r/6 Uia£j ; if j is labelled by 2 (12) 

e i4> (r/3 U 12 U 23 CJ , if j is labelled by 3 

on the original Hubbard Hamiltonian Hi 0 results in 
the new Hubbard Hamiltonian H 2 0 that shows the 
global U(l) spin-rotational symmetry clear ly.^1 Fig. i'i) 
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and (c) show the DM vectors in the original Hamlitonian 
Hi ([lj and the gauge-transformed Hamiltonian H 2 Q 
(with n = z), respectively. 


IV. GROUNDSTATES OF THE SOC SPIN 
HAMILTONIAN 

Starting from the Hubbard Hamiltonian _ff 2 |l]), the 
large U limit at half filling freezes the charge fluctua¬ 
tion, and eventually ends up with the spin Hamiltonian 
([ 5 ]) on the second order perturbation in t/t/PSEEI The 
SO(3) counterpart of the SU(2) gauge transformation in 
Eq. ([l2]), given by 

! S j , if j is labelled by 1 

R(— (/>n)i?.(di 2 )Sj, if j is labelled by 2 , 

R(— 20 n)i?(d 12 )i?(d 23 )Sj, if j is labelled by 3 

(13) 

yields the new spin Hamiltonian iff ©>■ The spin Hamil¬ 
tonian iff is invariant under global SO(2) spin rotation 
about the axis n, S( 1 — > R( 0 n) • S', which is the conse¬ 
quence of the U(l) spin-rotational symmetry of the par¬ 
ent Hubbard Hamiltonian. The full SO(3) spin-rotational 
symmetry is respected once the the flux vanishes <f> = 0 
(mod 2-7 t), or equivalently, 0 = 0, 27 t/ 3, or 47r/3PSl 
The dependence of the symmetry of the Hamiltonian 
on its parameter, 0 in our case, indicates the possible 
dramatic change in the physical properties when 0 cross¬ 
ing the high symmetry points 0 = 0, 27r/3, and 47 t/3 . To 
see that, we treat s pins classically in the spin Hamilto¬ 
nian fff ([ 5 ]) [27I47 48] which allows us to obtain the exact 
groundstates. The groundstates at the high symmetry 
points, where the fluxes are zero, are already known Jl2M9| 
We thus focus on nonvanishing fluxes. 

We start for an isolated triangle. Observing the Hamil¬ 
tonian is quadratic in spin, a straightforward way is to list 
all three spins into a large column spin S = (Si, S2, S3), 
and diagonalize a 9 x 9 matrix ff representing the spin 
Hamiltonian fff = S ■ ff ■ S/2. Specifically the matrix ff 
is given by 

/ 0 R R t \ 

U=\R t 0 R , (14) 

\R R t 0/ 

where R = f?( 0 n). This scheme is generally invalid as 
the eigenstates may not be physical due to the different 
lengths of spins, e.g., |Si| 7^ |S2 1- For the current prob¬ 
lem, however, the eigenvectors with the minimum energy 
turn out to satisfy |Si| = |S 2 | = | S3 1 always, which makes 
them physical. Spins are perpendicular to the SU(2) flux 
vector 4 >n in the continuously degenerate groundstates, 
which are given by 

f Sl = ff(^)S 2 = f?(Y)S3, o<0<f 
J Si = So = S 3 , f <0< f . (15) 

[Si = f?(2|h)S 2 =i?(4|h)s 3 , ^<0<2 tt 


For two neighboring triangles, their exact groundstates 
of the spin Hamiltonian can be patched by matching the 
spin of the shared site. This procedure can be extended 
to the entire Kagome lattice because all the triangles have 
the same DM vectors. Specifically, spins labeled by the 
same number [see Fig.jlja)] point in the same direction in 
the groundstates. For example, the groundstates of the 
spin Hamiltonian fff with 27r/3 < 0 < 47t/ 3 have all the 
spins pointing in the same direction in the plane perpen¬ 
dicular to n. Fig.[2](a) shows a groundstate of fff for the 
angle 0 < 0 < 2tt/3 and n = z. The exact groundstates 
of the spin Hamiltonian fff ([ 5 ]) can be obtained from 
those of fff (15) by reversing the gauge transformation 
in Eq. (13). 

The groundstates of the spin Hamiltonian break the 
SO(2) spin-rotational symmetry, which signals the exis¬ 
tence of the BKT transition at finite temperature. In the 
continuum approximation, the transition temperature is 
given by 


Tkt = 


(tta/ 3/4) S’ 2 cos(l>/3), (16) 


where 4> = $ mod 27 t, 2 t r < $ < 47 t, when neglect¬ 
ing spin waves. Fig. Hb) shows the results for the 
transition temperatures from finite temperature classi¬ 
cal Monte Carlo simulation of the spin Hamiltonian fff 
and various 0 with the spin length of unity, which agrees 
well with the theoretical prediction (16) with a renor¬ 


malized spin length S = 0.69. The renormalization of 
the spin length can be attributed to thermal spin-wave 
fluctuations. 


V. DISCUSSION 

We have showed that the single-band Hubbard Hamil¬ 
tonian on the Kagome has the global U(l) spin-rotational 
symmetry even in the presence of SOC. The U(l) sym¬ 
metry does not demand a specific shape of constituent 
triangles, but only requires the inversion symmetry be¬ 
tween neighboring triangles centered at the shared site. 
Linear deformation of the Kagome lattice by strain, for 
example, would break the three-fold rotational symme¬ 
try about the centers of the triangles, but preserves the 
inversion symmetry centered at sites and thus maintains 
the associated U(l) symmetry as well®! 

We have provided the exact classical groundstates of 
the spin Hamiltonian, which spontaneously break the 
continuous symmetry. The BKT transition occurs at fi¬ 
nite temperature, which is governed by the SU(2) flux 
threading the triangles. This is an example showing the 
physical effects of the SU(2) flux associated with SOC. 
Its effect on the quantum Hamiltonian would deserve to 
be investigated. We would like to mention that the con- 
trollabl e SU (2) gauge field has been created in optical 
lattices ElU2i which may afford the platform to observe 
the effects of varying flux on the physical properties of 
the quantum Hamiltonian. 
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Preservation of U(l) spin rotational symmetry depends 
on lattice structure. The Kagome lattice (with inversion 
symmetry centered at sites) is one of the lattices whose 
structure support U(l) spin rotational symmetry even 
in the presence of arbitrary SOC. The approach taken 
in this paper, however, can be applied to other lattice 
systems. For example, on the square lattice, the Hub¬ 
bard Hamiltonian with arbitrary SOC also possesses the 
global U(l) spin-rotational symmetry provided the inver¬ 
sion symmetry centered at sites is respected. The classi¬ 
cal groundstates of the corresponding spin Hamiltonian 
are the Neel states polarized along the SU(2) flux vector, 
for the square lattice is bipartite. The groundstates do 
not break the continuous symmetry, and thus the BKT 
transition does not occur. It would be worth pursuing 


to study other two- and three-dimensional lattices in the 
similar approach. 
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